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Electron- hole pairs produced by tunneling in a degenerate electron gas lose their spin entanglement 
by spin-orbit coupling, which transforms the fully entangled Bell state into a partially entangled 
mixed density matrix of the electron and hole spins. We calculate the dependence of the entangle- 
ment (quantified by the concurrence) on the spin-orbit coupling time r so and on the diffusion time 
(or dwell time) Tdweii of electron and hole in the conductors (with conductances ^> e 2 /h) at the 
two sides of the tunnel barrier (with conductance <C e 2 /h). The entanglement disappears when the 
ratio Tdwcii/Tso exceeds a critical value of order unity. The results depend on the type of conductor 
(disordered wire or chaotic quantum dot), but they are independent of other microscopic parameters 
(number of channels, level spacing). Our analytical treatment relies on an "isotropy approximation" 
(no preferential basis in spin space), which allows us to express the concurrence entirely in terms 
of spin correlators. We test this approximation for the case of chaotic dynamics with a computer 
simulation (using the spin kicked rotator) and find good agreement. 

PACS numbers: 71.70.Ej, 03.67.Mn, 72.25.Rb, 73.23.-b 



I. INTRODUCTION 

Spin-orbit coupling is one of the sources of degradation 
of spin entanglement that has been extensively investi- 
gated for electron pairs confined to two quantum dots. 1 
In that context the spin-orbit coupling induces dephas- 
ing by coupling the electron spins via the orbital motion 
to fluctuating electric fields in the environment (due to 
lattice vibrations or gate voltage fluctuations) . The cou- 
pling of the spins to the environment is needed for entan- 
glement degradation because the spin-orbit coupling by 
itself amounts to a local unitary transformation of the 
electron states in the two quantum dots, which cannot 
change the degree of entanglement. 

The characteristic feature of these quantum dots is 
that they are single-channel conductors with a conduc- 
tance G that is small compared to the conductance quan- 
tum e 2 /h. This implies in particular that the width of 
the energy levels is much smaller than the mean level 
spacing. At low voltages and temperatures there is then 
only a single accessible orbital mode. This is the main 
reason that spin-orbit coupling by itself cannot degrade 
the spin entanglement. 

In a multi-channel conductor the situation is altogether 
different. Flying qubits in a multi-channel conductor can 
lose their entanglement as a result of spin-orbit coupling 
even in the absence of electric field fluctuations, because 
the large number of orbital degrees of freedom can play 
the role of an environment. This mechanism is the elec- 
tronic analog of the loss of polarization entanglement by 
polarization-dependent scattering in quantum optics^^ 
Fully-phase-coherent spin-orbit coupling can degrade the 
spin entanglement by reducing the pure spin state to a 
mixed spin density matrix — which typically has less en- 
tanglement than the pure state. Here we investigate this 
mechanism in the context of electron-hole entanglement 
in the Fermi sea£ Apart from the practical significance 
for the observability of the entanglement, this study pro- 
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FIG. 1: Multi-channel conductor containing a tunnel barrier. 
The applied voltage creates electron-hole pairs (solid and open 
circles) at opposite side of the barrier, whose spin state is max- 
imally entangled. As the pair moves through the leads, the 
spin and orbital degrees of freedom become entangled by the 
spin-orbit coupling, degrading the spin entanglement upon 
tracing out the orbital degrees of freedom. 



vides a test for a theory of entanglement transfer based 
on the "isotropy approximation" that the spin state has 
no preferential quantization axis. 

The system we consider, a multi-channel conductor 
containing a tunnel barrier, is schematically depicted in 
Fig. [U The applied voltage V creates, at each tunnel 
event, a maximally entangled electron-hole pair£ Spin- 
orbit coupling in the leads entangles the spin and orbital 
degrees of freedom. The spin state (obtained by tracing 
out the orbital degrees of freedom) is degraded from a 
pure state to a mixed state. The degree of entanglement 
of the spin state decreases and can vanish for strong spin- 
orbit coupling. We consider two cases. In the first case 
the leads are diffusive wires while in the second case we 
model the leads as two chaotic cavities. Although the 
first case is our primary interest, we include the second 
case in order to test our approximate analytical calcula- 
tions against an exact numerical simulation of the spin 
kicked rotator^ 

The outline of the paper is as follows. In Sec. [IT] we 
calculate the density matrix of the electron-hole pairs 
in the regime where the tunnel conductance GtunneJ is 
•C e 2 /h. This is the regime that the electron-hole pairs 
form well separated current pulses, so that their entangle- 
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ment can be measured easily^ (For G tU imci ^ e 2 /h differ- 
ent eiectron-hoie pairs overiap in time, complicating the 
detection of the entanglement.) From the density matrix 
we seek, in Sec. IIII1 the degree of entanglement as mea- 
sured by the concurrence £ For our analytical treatment 
we approximate the density matrix by the spin-isotropic 
Werner stated The absence of a preferential basis in spin 
space is a natural assumption for a disordered or chaotic 
system, but it needs to be tested. For that purpose we 
use the spin kicked rotator, which is a stroboscopic model 
of a chaotic cavity^ We conclude in Sec. [TV] 



II. 



CALCULATION OF THE ELECTRON-HOLE 
STATE 

A. Incoming and outgoing states 



Since the scattering of both orbital and spin degrees 
of freedom is elastic, we may consider separately each 
energy E in the range (Ep , Ep + eV). For ease of notation 
we will omit the energy arguments in what follows. We 
assume zero temperature, so the incoming state is 
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where |0f) is the unperturbed Fermi sea, 
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|0 F ) = det(r) JJ b\ v |0> 



(5) 



Since rr^ — 1 — 0(e), we may assume that r is a unitary 
matrix to the order in e considered. The determinant 
det(r) is therefore simply a phase. The state ([!]) is a 
superposition of the unperturbed Fermi sea and a sin- 
gle electron-hole excitation, consisting of an electron in 
channel v at the right and a hole in channel /i at the left. 

As a check, we verify that the multi-channel result (H} 
reduces for N — 1 to the single-channel result 
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of Ref. |g. We use the identity^ 

a y r T a y — det(r)r\ 



(6) 



(7) 



which holds for any 2x2 unitary matrix r (with a y a Pauli 
matrix). Hence ta y r T a y = det(r)trt+C(e). Substitution 
into the single-channel result ((6|) indeed gives the multi- 
channel result ||3} for N = 1. 



The creation operators a L v , v = 1, . . . , 2N (acting on 
the true vacuum |0)) occupy the z^-th channel incoming 
from the left. The index v labels both the N orbital 
and two spin degrees of freedom. The 2N channels in- 
coming from the right (creation operators a R v ) are un- 
occupied in the energy range {E-p, Ep + eV). We collect 
the creation and annihilation operators in vectors &l = 

(a>L,i,a>L,2i ■ ■ ■ i cll,2n), &i? = (or,i, or,2, ■ ■ ■ , o,r,2n)- 

The annihilation operators bL, v and bn. v of the outgo- 
ing channels are related to those of the incoming channels 
by the scattering matrix 
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The AN x AN unitary scattering matrix S is decomposed 
into 27V x 2N transmission and reflection matrices t, t' , r, 
and r'. Substitution into Eq. ([1]) gives the outgoing state 



l*n,.t) = 
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B. Tunneling regime 

We expand the outgoing state ([3]) in the small param- 
eter e = (/i/e 2 )Gtunnci, neglecting terms of order e and 
higher. Since t,t' are 0(e 1 ^ 2 ) while r,r' are C(e°), we 
keep only terms linear in t and t' . The result is 

l^out) = |0p) + E^ rt WOi.M IOf) + O(e), (4) 



C. Spin state of the electron- hole pair 

The spin state of the electron-hole pair is obtained from 
l^out) by projecting out the vacuum contribution and 
then tracing out the orbital degrees of freedom. This 
results in the 4x4 density matrix 
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with w = tr (tHr'r). Here n and m label the orbital de- 
grees of freedom and a, (3, 7, and 5 label the spin degrees 
of freedom. 

We assume that the tunnel resistance is much larger 
than the resistance of the conductors at the left and right 
of the tunnel barrier. The transmission eigenvalues T n 
(eigenvalues of W) are then determined mainly by the 
tunnel barrier and will depend only weakly on the mode 
index n. We neglect this dependence entirely, so that 
T n = T for all n, the tunneling conductance being given 
by Gtunnd = {2e 2 /h)NT. 

To obtain a simpler form for the density matrix we use 
the polar decomposition of the scattering matrix 



S = 



r t' 
t r' 



u 
v 



u' 
v' 



(9) 

where u, v! , v, and v' are unitary matrices and T = 
diag(Ti, T2, . . . , T2n)- For mode independent T„'s the 
matrix T equals T times the unit matrix. Hence 



tr 1 = y/(l-T)TU 



(10) 
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is proportional to the 2N x 27V unitary matrix U = vvf. 
Substitution into the expression ([5]) for the density ma- 
trix gives 



1 

2N 



Una,mf3U nlmS . 



(11) 



If there is no spin-orbit coupling, the matrix U is di- 
agonal in the spin indices: U na ,mSi — U nm 8 a p with U an 
N x N unitary matrix. The density matrix then repre- 
sents the maximally entangled Bell state IV'Bell)) 

(pBell)a,8,7<5 = r^a/S&y* = I^Bcll) (V'BollI , (12) 
IV-Bcll) =^(|T)elT)h + l4) e li)h)» (13) 

with \cr) h an electron (e) or hole (h) spin pointing up 
(a =1") or down (a =],). The state (fP2"]) is a pure state 
(PBeil = PBeil)' Spin-orbit coupling will in general de- 
grade p to a mixed state, with less entanglement. 



III. ENTANGLEMENT OF THE 
ELECTRON-HOLE PAIR 

We quantify the degree of entanglement of the mixed 
electron-hole state (j8]) by means of the concurrence C 
(which is in one-to-one correspondence with the entangle- 
ment of formation and varies from for a nonentangled 
state to 1 for a maximally entangled state). Following 
Wootters^ the concurrence is given by 



C = max 



{0, 



Ai - v/Aa - v As — v A 



} , (14) 



where the Ai's are the eigenvalues, in decreasing order, 
of the matrix product p(a v (g> a y )p*{a y ® a y ). 

In the next two subsections we calculate the concur- 
rence numerically for a chaotic cavity using Eq. {§]) and 
analytically with an isotropy approximation for the den- 
sity matrix. 



A. Numerical simulation 

We calculate the concurrence C numerically for the 
case that the scattering at the left and at the right of 
the tunnel barrier is chaotic. (The more experimentally 
relevant case of diffusive scattering will be considered in 
the next subsection.) 

The total scattering matrix S of the system (shown in 
Fig. is constructed from the scattering matrix of the 
tunnel barrier, 



St 



N 



N 



(15) 



and the scattering matrices Si and 5*2 of the cavity on 
each side of the tunnel barrier. (We denote by Ijv the 
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FIG. 2: Two chaotic cavities with scattering matrices Si and 
52 connected by a tunnel barrier with scattering matrix St- 
The chaotic cavities are modeled by two spin kicked rotators. 



N x N unit matrix.) We expand S in the small param- 
eter T and keep terms up to order ©(T 1 / 2 ) = C(e 1/2 ), 
consistent with the expansion of the outgoing state ([4]). 
This results in 
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r = r 1 +t' 1T —-jt 1 +0(T), 



t = t? 



1 



T- 



l + r 2 1 



(16a) 
(16b) 



and similar expressions for r' and t' which we do not 
need. 

The scattering matrices S\ and S% of the chaotic cav- 
ities are constructed from two spin kicked rotators. We 
briefly describe this method of simulation in Appendix [Al 
referring to Ref. H for a more detailed exposition. 

The resulting ensemble-averaged concurrence as a 
function of the ratio Td W0 n/T so of the mean dwell time 
Tdwcii and spin-orbit coupling time r so is shown in Fig. [3] 
The dwell time Tdweii is the average time between a tunnel 
event and the escape of the particle into the left or right 
reservoir. The time t so is the exponential relaxation time 
of the spin-up and spin-down densities towards the equi- 
librium distribution. (Both time scales are calculated in 
Appendixes [Al and IBP For a single channel, N = 1, the 
concurrence is unity independent of spin-orbit coupling 
strength since the trace over the orbital degrees of free- 
dom leaves p unchanged. From Fig. [3] (bottom panel) 
we see that for small N the concurrence saturates at a 
nonzero value for large Tdweii/ T so : 



lim 

ell/Tso- 



(C> = 



1, 


N = 


1, 


0.15, 


N = 


2. 


0.01, 


N = 


3. 



(17) 



The limiting value for N = 2 is close to that obtained 
in Ref. [l^ in a single chaotic cavity. For TV > 5 
the ensemble-averaged concurrence is negligible for large 
Tdweii/Tgo- The dependence of (C) on Td we n/T so becomes 
N independent for N > 15. 
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FIG. 3: Ensemble averaged concurrence of the electron-hole 
pair scattered by two chaotic cavities, as a function of the 
spin-orbit coupling rate l/r so for different number of modes, 
N, in the leads. The data points are calculated from the spin 
kicked rotator; the lines are guides to the eye. The middle 
panel shows that the results become iV-independent for large 
N while the bottom panel shows that for small N the concur- 
rence saturates at a finite value. 



B. Isotropy approximation 



(v ® v*)p(y* ® v T ) = p 



(18) 



for all 2 x 2 unitary matrices V. This transformation 
rotates the spin basis of the electron (acted on by V) 
and the hole (acted on by V*) by the same rotation angle. 
The isotropy relation (|18p constrains the density matrix 
to be of the Werner form 



PW = i(l-£) 1 4 + £ IV'Bell) (V>Bell| 



1 



< t < 1, (19) 



with IV'Bell) the Bell state defined in Eq. (fT3| . The con- 
currence of the electron-hole Werner state is given by 



C(p 



w ) 



3/1 
-maxICU-- 



(20) 



The parameter £ characterizing the Werner state can 
be calculated from 



£ = tr [{a z <g) a z )p] = p n - P22 - P33 + Pa- 



(21) 



Only diagonal elements of the density matrix appear in 
the expression (|2"Tj) for £. These can be calculated semi- 
classically in the TV-independent limit N ^> 1 (see Ap- 
pendix [BJ, leading to the following expressions for the 
concurrence: 



(C) 

diffusive 

(C) 

chaotic 



0, < T so < 1.5 ^dwcll 1 

47r(-l)"(2n+ 1) 



7r 2 (2n+ l) 2 +8r dwen /T s , 



(22a) 



r(l + Tdwcii/Tso) 2 - h 



0. 



'dwell <~ _ 

_>■ ^3-1 ^ so ' 

n < r < Tdwc " 



(22b) 



In Fig. [4] we plot the analytical result (|22|) for the con- 
currence. The two cases of diffusive and chaotic scatter- 
ing differ only slightly. The initial slopes are the same, 

(C) 

diffusive — 

(C) 

chaotic 

(23) 

= 1 - 3r dwc ii/T so + C(r dwcl i/T so ) 2 . 

The critical spin-orbit coupling strengths, beyond which 
the concurrence vanishes, are different: T™ tlcal = 
l-5rdwcii for diffusive scattering and r™ tlcal = 
T'dweii/(v / 3 — 1) = 1.37r,i we ii for chaotic scattering. 

We also compare in Fig. [4] the analytical results in the 
chaotic case from this section with the numerical results 
from the previous section. The agreement is quite good 
for large N, where the semiclassical analytics is expected 
to hold. 



To obtain an analytical expression for the entangle- 
ment degradation we approximate the density matrix by 
the spin-isotropic Werner state.— The absence of a pref- 
erential basis in spin space means that the density matrix 
p for an electron-hole pair is invariant under the trans- 



IV. CONCLUSION 

Figure U summarizes our main findings: The effect of 
spin-orbit coupling on the degree of spin-entanglement of 
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FIG. 4: Ensemble averaged concurrence as a function of spin- 
orbit coupling rate l/r so . The solid and dashed curves are the 
analytical results (1221! in the case of diffusive wires and chaotic 
cavities, respectively, on each side of the tunnel barrier. These 
analytical curves use the isotropy approximation. The data 
points are from the numerical simulation in the chaotic case 
without the isotropy approximation (spin kicked rotator of 
Fig. El with N = 30). 



so far not been accomplished experimentally. Thus, the 
isotropy approximation that has been used here as a way 
to simplify the calculation of the concurrence, also has an 
experimental implicationjlS, By relying on spin isotropy 
the concurrence can be obtained directly from correla- 
tors of time averaged spin currents. Our demonstration 
of the accuracy of the isotropy approximation may mo- 
tivate experimentalists to try this "poor man's method" 
of entanglement detection — since average spin currents 
have been measured^ while spin noise has not. 
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APPENDIX A: SPIN KICKED ROTATOR 



the electron-hole pairs produced at a tunnel barrier de- 
pends strongly on the ratio of the dwell time Td wo n and 
spin-orbit coupling time t so . Even though Td wc ii and r so 
each depend sensitively on the nature of the dynamics 
(diffusive or chaotic) the dependence of the concurrence 
on the ratio Tdweii/^so is insensitive to the nature of the 
dynamics. The initial decay (|23|) is the same and the 
critical spin-orbit coupling strength (beyond which the 
entanglement vanishes) differs by less than 10%. This 
has the useful experimental implication that a single pa- 
rameter suffices to quantify the amount of entanglement 
degradation by spin-orbit coupling. 

We have tested our analytical theory using a com- 
puter simulation for the case of chaotic dynamics. (The 
close similarity to the diffusive results suggests that this 
test is representative.) Analytics and numerics are in 
good agreement, differing by less than 10% in the regime 
N ^> 1 of large conductance G where the semiclassical 
analytics applies. While the semiclassical approximation 
is controlled by the small parameter 1/N (or, more gener- 
ally, h/e 2 G), the isotropy approximation has no small pa- 
rameter that controls the error. Its use is justified by the 
reasonable expectation that an ensemble of disordered or 
chaotic systems should have no preferential quantization 
axis for the electron and hole spins. It is gratifying to see 
that the numerics supports this expectation. 

The standard method of experimentally verifying the 
presence of entanglement is by demonstrating violation of 
Bell inequalities. In optics this is achieved by measuring 
coincidence rates of photons by photodetectors (i.e. by 
counting photons) in different polarization bases. In the 
solid state one cannot simply count electrons, but rather 
needs to formulate the Bell inequalities in terms of cor- 
relators of spin currents (= spin noise)i 13 ' 14 i 15 This has 



We use the spin kicked rotator to model transport in 
chaotic cavities in the presence of spin orbit coupling. 
The model is described in detail in Ref. ,8. Here we collect 
the most important results of that paper for reference. 

The spin kicked rotator is a numerically efficient stro- 
boscope description of the chaotic dynamics, described 
by a quantum map with Floquet matris*^ 



{JLUXU^U) W , 1,1' = 0,1,... ,M -1. (Al) 



The matrices appearing in this expression are defined by 

H„, = 5 w e-^ l+l ^l M l 2: (A2a) 
U w = A/- 1/2 e- j2 ^'/ M l 2 , (A2b) 

X W = 5we -i{M/A*)V{2*l/M)^ (A2c) 

with kicking potential 

V{6) = K cos6 t 2 + K so {a x sin26 + cr z sin8). (A3) 

The matrices a x and a z are Pauli matrices. The dimen- 
sion 2M x 2M of the Floquet matrix determines the 
mean spacing 8 = 2ir/M of the quasienergies e. The 
map is classically chaotic for kicking strength K > 7.5. 
The parameter K so breaks spin rotation symmetry and 
lo breaks other symmetries of the map. In our simula- 
tions we choose K = 41 (fully chaotic), M = 640, and 
Iq = 0.2. The spin-orbit coupling time r so (in units of the 
stroboscopic period) in the model is determined by the 
parameter A^ so through^ 



16tt 2 
KIM 2 ' 



(A4) 



The reason for denoting the variable of the kicking po- 
tential by 9 is historical: In its original incarnation the 
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FIG. 5: Location in classical phase space of the open bound- 
ary conditions representing the leads. 



Tdweii ( in units of the stroboscopic period) is given by 



M _ 2tt 

Tdwell - — - — , 



(A7) 



where we have taken into account the fact that N of the 
2N channels are closed by the tunnel barrier (cf. Fig. 
Notice that the mean dwell time is a classical quantity, 
while N and 5 separately are quantum mechanical quan- 
tities. 

From the reflection and transmission matrices (|I6p the 
density matrix © is obtained, from which the concur- 
rence (|14p follows. The concurrence is averaged over 20 
different quasienergies e, e' and over 20 different lead 
positions P, P' in the two cavities (assumed to be inde- 
pendent scatterers) . Results are shown in Fig. [3] 



model was a freely moving rotator (pendulum) periodi- 
cally kicked by the kicking potential V with K so = 0J£ 
The Floquet matrix was then derived from the Hamilto- 
nian of this kicked rotator. For our purposes it is much 
more fruitful to consider the Floquet matrix (time evolu- 
tion operator) T as the definition of the model. One then 
considers the matrix X to describe the free (spin-orbit 
coupled) motion inside the chaotic cavity. This free mo- 
tion is interrupted by boundary scattering which is given 
in terms of the matrix IL (The matrix X is diagonal in 9 
space, while II is diagonal in p-space; U maps between the 
two spaces. 9 and p are conjugate variables.) With this 
interpretation, the variable 9 becomes the momentum- 
like variable, while p becomes a variable for the position 
on the boundary. To accommodate the notation it can be 
useful to think of 9 as the angle describing the direction 
of the momentum. 

The Floquet matrix T represents the internal dynam- 
ics of the chaotic cavity. In order to use the map to 
describe an open cavity, we need to open up the map. 
This is achieved by imposing open boundary conditions 
in a subspace of Hilbert space; that is, any time the map 
takes the particle to this subspace it is absorbed. With 
the above interpretation in mind, it becomes intuitively 
clear that these open boundary conditions take the form 
of strips parallel to the 9 axis in phase space (cf. Fig. [5} . 
We thus represent the boundary conditions by the indices 
for the first lead and for the second lead (the sub- 
script k — 1,2, ... ,N labels the modes). These indices 
label the momentum eigenstates defining the subspace. 
The 4N x 2M projection matrix 



Pka,k'f3 



5 a 8 if k' = I 



i _ j(aO 







otherwise, 



(A5) 



projects onto the leads. The scattering matrices of the 
two chaotic cavities now take the for m 20 i 21 i 22 i 23 

S 1 =P[e~ iE -F(l-P T P)]- 1 J r P T , (A6a) 
S 2 = P'[e- lE ' - T{\ - P' T P')]- 1 TP' T , (A6b) 

where P and e refer to the first chaotic cavity and P' and 
e' refer to the second chaotic cavity. The mean dwell time 



APPENDIX B: CALCULATION OF SPIN 
CORRELATORS 

The diagonal elements of the density matrix appear- 
ing in the expression (|2"Tj) for the Werner parameter £ 
represent spin correlators, 

Pn = P^h P22 = P]li P33 = p lh P44 = Pa- (Bl) 

Here P aa i is the probability that the outgoing electron 
has spin a and the outgoing hole has spin a'. To calculate 
these correlators, it is convenient to first consider only 
those electrons that exit after a time t and those holes 
that exit after a time t'. The time- resolved correlator 
P aa '(t,t') gives the desired P aa * after integration over 
time, 

/>oo />oo 

P**>= dt dt'P aa ,(t,t')P dwcll {t)P dwcU (t'), (B2) 
Jo Jo 

weighted by the dwell time distribution Pdweii (we assume 
that the dwell times at the left and right of the tunnel 
barrier are independent and identically distributed). 
As initial condition we take 



p a<r .(o,o) = -<W, 



(B3) 



corresponding to the spin state immediately after the 
tunnel event. Spin-orbit coupling randomizes the spin 
with a rate l/r so , so that P aa '{t,t') decreases in time 
according to the rate equations 

^-P^(t,t') =-L^[P CT , v ,(t,t')-P ffCT > (t,t% (B4a) 



dt 

d 
di 



7 P**>(t,t') 



' ^[P CTCT „(M') -P aa , {t,t')} . (B4b) 



2t so 



The solution of the rate equations (|B4|) with the initial 
condition (|B3|) is 

P TT (M') = Pu(M') = \ + \e-^/^, (B5a) 
P u (t,t') = P u {t,t') = \- I e -('+*')M„. (B5b) 
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FIG. 6: Dwell time distribution in a diffusive wire (solid line) 
and chaotic cavity (dashed line). 



To complete the calculation we need the dwell time 
distribution. For a chaotic cavity this has the well known 
exponential forr^' 



P, 



1 



dwell, chaotii 



T dwell 



-</ r dwcll 



with 



Tdwell 



2ttH 
iVA 



(B6) 



(B7) 



inversely proportional to the mean level spacing A of 
Kramers degenerate levels in the cavity. 

For the diffusive wire (diffusion constant D) we de- 
termine Pdweii by solving the one-dimensional diffusion 
equation 



with initial and boundary conditions 

^(0,t)=0, p(L,t) = 0, p(x,0) = S(x). (B9) 

Here p(x,t) is the classical probability of finding a par- 
ticle at point x at time t. The boundary conditions rep- 
resent reflection by the high tunnel barrier at x = and 
absorption by the reservoir at x = L. 

The probability that the particle is still in the wire at 
time t is given by 



N(t) = / p(x,t)dx, 
Jo 

and therefore the dwell time distribution is 

dN(t) 



-Pdwcll — — " 



dt 



(BIO) 



(Bll) 



Solution of the diffusion equation by expansion in eigen- 
states gives the result in the form 



-Pdwell, diffusive = 7; > (— l) n (2n+l 

^Tdwell ' 



-(2n+l) 2 



71=0 



The mean dwell time is 



T dwell 



2D' 



(BI2) 



(B13) 



The dwell time distributions for the chaotic and diffusive 
dynamics are compared in Fig. [6l 

Collecting results we arrive at the expressions (f2"2"j) for 
the concurrence given in the main text. 
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